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Abstract 


The normalized differential cross section for top quark pair (tt) production is mea¬ 
sured in pp collisions at a centre-of-mass energy of 8 TeV at the CERN LHC using 
the CMS detector in data corresponding to an integrated luminosity of 19.7fb^^. The 
measurements are performed in the lepton+jets (e/ fi+iets) and in the dilepton (e+e^, 
, and decay channels. The tt cross section is measured as a function of 

the kinematic properties of the charged leptons, the jets associated to b quarks, the 
top quarks, and the tf system. The data are compared with several predictions from 
perturbative QCD up to approximate next-to-next-to-leading-order precision. No sig¬ 
nificant deviations are observed relative to the standard model predictions. 
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1 Introduction 

Understanding the production and properties of top quarks is fundamental for testing the qual¬ 
ity of the standard model (SM) and for searching for new physical phenomena beyond its scope. 
The large top quark data samples produced in proton-proton (pp) collisions at the CERN LHC 
provide access to precision measurements that are crucial for checking the internal consistency 
of the SM at the LHC energy scale. In particular, measurements of the top quark pair (tf) 
production cross section as a function of tf kinematic observables are important for compar¬ 
ing with the state-of-the-art quantum chromodynamic (QCD) predictions within the SM, and 
thereby constrain QCD parameters. In addition, the top quark plays a relevant role in theories 
beyond the SM, and such differential measurements are therefore expected to be sensitive to 
new phenomena (ll. 

Differential tf production cross sections have been measured previously at the Lermilab pp 
Tevatron [SO, and at the LHC at a centre-of-mass energy = 7TeV HHSI. We present 
here the first measurement of the normalized differential tf production cross section with the 
CMS detector at v^s = 8TeV. The analysis uses data recorded in 2012 corresponding to an 
integrated luminosity of 19.7 ± 0.5 fb^^, which is about a factor of four larger than the sample 
used in the measurement performed by the CMS Collaboration at 7TeV |5l. The analysis largely 
follows the procedures of Ref. [O and benefits from the increase in statistical precision together 
with improvements in kinematic reconstruction algorithms and extended systematic studies, 
leading to a significant reduction of the total uncertainties. 

The measurements are performed in f+jets channels (f = e or ji), which contain a single iso¬ 
lated charged lepton and at least four jets in the final state, and in dilepton channels, with two 
oppositely charged leptons (e+e^, and at least two jets. The tf cross section is 

determined as a function of the kinematic properties of the top quarks and of the tf system, as 
well as of the leptons and jets associated with bottom (b) quarks (b jets) from top quark decays. 

The kinematic properties of top quarks are obtained through kinematic-fitting and reconstruc¬ 
tion algorithms. The normalized differential tf cross section is determined by counting the 
number of tf signal events in each bin of a given observable, correcting for detector effects and 
acceptance, and dividing by the measured total inclusive tf event rate. The latter is evaluated 
by integrating over all bins in each observable. 

The results for directly measured quantities, such as kinematic properties of leptons and b jets, 
are presented in a fiducial phase space defined by the kinematic and geometric acceptance of 
all selected final-state objects. This avoids extrapolating the measured cross section into re¬ 
gions that are not experimentally accessible. In addition, the top quark and tf distributions are 
determined in the full phase space, in order to facilitate the comparison with higher-order per¬ 
turbative QCD calculations. The results are compared to several predictions obtained with the 
leading-order (LO) MadGraph |7] generator interfaced to pythia |8l for parton evolution and 
hadronization, the next-to-leading-order (NLO) generators powheg EHTTl, interfaced to both 
PYTHIA and HERWIG |[l2l , and MC@NLO llT^ interfaced to HERWIG, and the latest NLO calcu¬ 
lations with next-to-next-to-leading-logarithm (NNLL) corrections Ifl4l HSlI . and approximate 
next-to-next-to-leading-order (NNLO) predictions IIT^ . The approximate NNLO predictions 
can be computed with the DifeTop ITT I program. 

This document is structured as follows. A brief description of the CMS detector is provided 
in Section Details of the event simulation are given in Section and event reconstruction 
and selection are discussed in Section The estimated systematic uncertainties on the mea¬ 
surements of the cross section are described in Section The results of the measurement are 
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3 Event simulation and theoretical calculations 


discussed in Section]^ followed by a summary in Section]^ 

2 CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 13 m length and 6 m 
inner diameter, which provides an axial magnetic field of 3.8 T. Within the field volume are 
a silicon-pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), 
and a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two 
endcap sections. Charged particle trajectories are measured by the inner tracking system, cov¬ 
ering a pseudorapidity range of \i]\ < 2.5. The ECAL and the HCAL surround the tracking 
volume, providing high-resolution energy and direction measurements of electrons, photons, 
and hadronic jets up to |//| <3. Muons are measured in gas-ionization detectors embedded 
in the steel flux return yoke outside the solenoid covering the region \i]\ < 2.4. Extensive for¬ 
ward calorimetry complements the coverage provided by the barrel and endcap detectors up 
to \tj\ < 5.2. The detector is nearly hermetic, allowing for energy balance measurements in the 
plane transverse to the beam directions. A two-tier trigger system selects the pp collisions for 
use in the analysis. A more detailed description of the CMS detector, together with a definition 
of the coordinate system and the relevant kinematic variables, can be found in Ref. ITSl . 


3 Event simulation and theoretical calculations 

Event generators, interfaced with a detailed detector simulation, are used to model experi¬ 
mental effects, such as consequences of event reconstruction and choice of selection criteria, 
as well as detector resolution. The tt sample is simulated using the LO MadGraph event 
generator (v. 5.1.5.11), which implements the relevant matrix elements with up to three addi¬ 
tional partons. The MadSpin IIT^ package is used to incorporate spin correlation effects with 
matrix elements for up to three additional partons. The value of the top quark mass is fixed 
to nit = 172.5 GeV and the proton structure is described by the parton distribution functions 
(PDE) GTEQ6L1 Il20l . The generated events are subsequently processed with PYTHIA (v. 6.426, 
referred to as PYTHIA6 in the following) for parton showering and hadronization, and the MEM 
prescription Il21]| is used for matching of matrix-element jets to parton showers. The GMS de¬ 
tector response is simulated using Geant4 (v. 9.4) Il22ll . 

In addition to the MadGraph prediction, calculations obtained with the NLO generators 
MC@NLO (v. 3.41) and POWHEG (v. 1.0 rl380) are compared to the results presented in Sec¬ 
tion 1^ While POWHEG and MC@NLO are formally equivalent up to the NLO accuracy, they 
differ in the techniques used to avoid double counting of radiative corrections that can arise 
from interfacing with the parton showering generators. Two POWHEG samples are used: one is 
processed through PYTHIA6 and the other through HERWIG (v. 6.520, referred to as HERWIG6 in 
the following) for the subsequent parton showering and hadronization. The parton showering 
in PYTHIA6 is based on a transverse-momentum-ordered evolution scale, whereas in HERWIG6 
it is angular-ordered. The events generated with MC@NLO are interfaced with HERWIG6. The 
HERWIG6 AUET2 tune 1231 is used to model the underlying event in the POWHEG+HERWIG6 
sample, while the default tune is used in the mc@NLO+herwig6 sample. The proton structure 
is described by the PDE sets GTIO Il24ll and GTEQ6M 1201 for POWHEG and MC@NLO, respec¬ 
tively. In addition, the latest available NLO+NNLL IIT4l ITSl and approximate NNLO QGD 
predictions IT^ are also used to compare with the data. The NNLO MSTW2008 Il25l PDE set is 
used for both the NLO+NNLL and the approximate NNLO calculations. 

Standard model background samples are simulated with MadGraph (without the MadSpin 
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package), powheg, or pythia6, depending on the process. The main background contribu¬ 
tions originate from the production of W and Z/ 7 * bosons with additional jets (referred to as 
W+jets and Z+jets, respectively, in the following), single top quark (s-, t-, and tW channels), 
diboson (WW, WZ, and ZZ), tf production in association with a Z, W, or 7 boson (referred to 
as tt +Z/W /7 in the following), and QCD multijet events. The W+jets, Z+jets, and tt +Z/W /7 
samples are simulated with MadGraph with up to two additional partons in the final state. 
The POWHEG generator is used for simulating single top quark production, while PYTHIA6 is 
used to simulate diboson and QCD multijet events. Parton showering and hadronization are 
also simulated with PYTHIA 6 in all the background samples. The PYTHIA 6 Z2* tune 1261 is used 
to characterize the underlying event in both the tt and the background samples. 

For comparison with the measured distributions, the event yields in the simulated samples are 
normalized to an integrated luminosity of 19.7fb^^, according to their predicted cross sections. 
These are taken from NNLO (W+jets 1271128| and Z+jets ETI T NLO+NNLL (single top quark 
S-, t-, and tW channels fl^ T NLO (diboson 1291 . tt +W [30|, and tt +Z I13T1I 1. and LO (QCD mul¬ 
tijet 0) calculations. The predicted cross section for the tf +7 sample is obtained by scaling the 
LO cross section obtained with the Whizard event generator 1321 by an NLO/LO correction 
fC-factor 1331 . Correction factors described in Sections]^ andand subsequently referred to as 
scale factors, are applied when needed to improve the description of the data by the simula¬ 
tion. The tf simulation is normalized to the data to present the expected rates in the figures in 
Section m 

4 Event reconstruction and selection 

The event selection is similar to that described in Ref. |5j for the measurement of normalized 
differential tf cross sections at ^ys = 7 TeV, and is based on the final-state topology of tf events. 
The top quark decays almost exclusively into a W boson and a b quark, and only the subsequent 
decays of one or two of the W bosons into a charged lepton (electron or muon) and a neutrino 
are considered. These signatures imply the presence of isolated leptons with high transverse 
momentum pj, large pj imbalance caused by the neutrinos that escape detection, and highly 
energetic jets. The identification of b jets through b-tagging techniques is used to increase the 
purity of the selected sample. The event selection in each channel is optimized to maximize the 
content of tf signal events and background rejection. 

4.1 Lepton, jet, and missing transverse energy reconstruction 

Events are reconstructed using a particle-flow technique Il34ll35l . which combines signals from 
all subdetectors to enhance the reconstruction and identification of individual particles ob¬ 
served in pp collisions. Charged hadrons from pileup events, i.e. those originating from ad¬ 
ditional pp interactions within the same bunch crossing, are subtracted on an event-by-event 
basis. Subsequently, the remaining neutral-hadron component from pileup is accounted for 
through jet energy corrections 1361 . 

Electron candidates are reconstructed from a combination of the track momentum at the main 
interaction vertex, the corresponding energy deposition in the ECAL, and the energy sum of 
all bremsstrahlung photons attached to the track 1371 . The candidates are required to have 
Pj > 33GeV within the pseudorapidity interval \r]\ <2.1 for the £+jets channels, while elec¬ 
tron candidates in the dilepton channels are required to have pj > 20GeV and \ri\ < 2.4. 
As an additional quality criterion, a relative isolation lrei(0-3) < 0.10 in the ^+jets channels 
and lrei(0.3) < 0.15 in the dilepton channels is required, where Irei(^) is defined as the sum 
of the Pt of all neutral and charged reconstructed particle candidates inside a cone of AR = 
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4 Event reconstruction and selection 


'J + (A^)^ < X around the electron (excluding the electron itself) in rj-cp space, divided 
by the pj of the electron. 

Muon candidates are reconstructed using the track information from the silicon tracker and 
the muon system. They are required to have pj > 33GeV and \t]\ < 2.1 in the f+jets channels, 
while in the dilepton channels the corresponding selection requires pj > 20 GeV and \p\ < 2.4. 
Isolated muon candidates are selected if they fulfill Jj.ei(0.4) < 0.12 and 4ei(0.3) < 0.15 in the 
£+jets and dilepton channels, respectively. The same definition of relative isolation described 
above is also used for muon candidates. 

Jets are reconstructed by clustering the particle-flow candidates 1381 using the anti-kj cluster¬ 
ing algorithm with a distance parameter of R = 0.5 Il39l . Electrons and muons passing less 
stringent selections on lepton kinematic quantities and isolation, relative to the ones specified 
above, are identified but excluded from clustering. A jet is selected if it has pj > 30 GeV 
and \r]\ < 2.4 for both the ^+jets and dilepton channels. Jets originating from b quarks are 
identified through a "combined secondary vertex" algorithm ||40| , which provides a b-tagging 
discriminant by combining secondary vertices and track-based lifetime information. The cho¬ 
sen working point in the ^+jets channels has an efficiency for tagging a b jet of ~60%, while the 
probability to misidentify light-flavour jets as b jets (mistag rate) is only «1.5%. In the dilepton 
channels, the working point is selected to provide b-tagging efficiency and mistag rate of ~80- 
85% and ~10%, respectively HOI . These requirements are chosen to reduce the background 
contribution in the corresponding channels while keeping a large fraction of the tf signal. 

The missing transverse energy Ex is defined as the magnitude of the imbalance in the transverse 
momentum in the event, which is the negative of the vectorial sum of the momenta in the 
transverse plane of all the particles reconstructed with the particle-flow algorithm tffl . To 
mitigate the effect of contributions from pileup on the resolution in Ex/ we use a multivariate 
correction where the input is separated into components that originate from the primary and 
other collision vertices |42|. This correction improves the Et resolution by ft!5%. 


4.2 Event selection 


Events in the £+jets channels that are triggered by the presence of a single electron (muon) with 
Pj > 27GeV (pj > 24GeV, \rj\ < 2.1), are selected if they contain exactly one reconstructed 
lepton fulfilling the requirements described in Section [4T| Events are rejected if there are addi¬ 
tional electron candidates with pj > 20 GeV, \rj\ < 2.5, and Jrei(0.3) < 0.15, or additional muon 
candidates with pj > 10 GeV, \rj\ < 2.5, and J].ei(0.4) < 0.2. Additionally, an event must con¬ 
tain at least four reconstructed jets satisfying the criteria described in Section 4.1 To suppress 
background contribution mainly from W+jets events, at least two of these jets are required to 
be tagged as b jets, and at least two must not be tagged as b jets, as they are used to reconstruct 
W —> qq' decays. In the dilepton channels, events are triggered using combinations of two 
leptons with pj thresholds of 8 and 17 GeV, and are selected if they contain at least two isolated 
leptons of opposite electric charge and at least two jets. At least one of the jets is required to 
be b-tagged. In events with more than two leptons, we choose the lepton pair with opposite 
charge and largest value in the sum of their scalar pj. Events with an invariant mass of the 
lepton pair smaller than 20 GeV are removed to suppress events from decays of heavy-flavour 
resonances and low-mass Drell-Yan processes. Backgrounds from Z+jets processes in the e+e^ 
and channels are also suppressed by requiring the dilepton invariant mass to be outside 
a Z boson mass window of 91 ± 15 GeV, and to have Et > 40 GeV. 


After these selection steps, several basic distributions in £+jets and dilepton events are shown 
in Eigs.[^and|^ respectively. The hatched regions correspond to the shape uncertainties for the 
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signal and background (cf. Section and are dominated by the former. The data are reason¬ 
ably well described by the simulation, as shown in the lower part of each plot, where the ratio 
of data to simulation is presented to better indicate the level of agreement between data and 
the default tf signal (MadGraph+pythia6) and background samples used in the analysis. For 
both charmels, however, data tend to have lower pj values than predicted by the simulation. 
It has been verified that the results presented in Section are not affected by these remaining 
differences between data and simulation. A better data-to-simulation agreement in the lepton 
and jet pj distributions is obtained by scaling the top quark pj spectrum in simulation to match 
the data. However, the impact on the measurement of the cross sections is negligible. 
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Figure 1: Kinematic distributions after event selection and before the kinematic reconstruction 
of the tf system in the ^+jets channels: the multiplicity in the reconstructed number of b-tagged 
jets (top left), the multiplicity in the reconstructed number of jets (top right), the pj of the 
selected isolated leptons (bottom left), and the pj of all reconstructed jets (bottom right). The 
QCD multijet background is negligible and not shown. The hatched regions correspond to the 
shape uncertainties for the signal and backgrounds (cf. Section |^. The lower part of each plot 
shows the ratio of data to the predictions. 
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4 Event reconstruction and selection 
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Figure 2: Kinematic distributions after event selection and before the kinematic reconstruction 
of the tt system for the dilepton channels: the multiplicity in the reconstructed number of b- 
tagged jets (top left), the multiplicity in the number of reconstructed jets (top right), the pj of 
the selected isolated leptons (bottom left), and the px of the reconstructed jets (bottom right). 
The QCD multijet background is negligible and not shown. The Z/ 7 *+jets background is de¬ 
termined from data ||5H43l. The hatched regions correspond to the shape uncertainties for the 
signal and backgrounds (cf. Section]^. The lower part of each plot shows the ratio of data to 
the predictions. 


























































4.3 Kinematic reconstruction of the tt system 


7 


4.3 Kinematic reconstruction of the tt system 

The kinematic properties of the top quark pair are determined from the four-momenta of 
all final-state objects through kinematic reconstruction algorithms. These algorithms are im¬ 
proved versions of those described in Ref. (S]. 

In the £+jets channels, a constrained kinematic fitting algorithm is applied OIHl to the four- 
momenta of the selected lepton and up to five leading jets, and the representing the trans¬ 
verse momentum of the neutrino, which are changed according to their resolutions. The fit is 
constrained to reconstruct two W bosons, each with a mass of 80.4 GeV. In addition, the recon¬ 
structed top quark and antiquark masses are required to be equal. To reduce the number of 
permutations in the association of jets to quarks, only b-tagged jets are considered as b quarks, 
and only untagged jets are considered as light quarks. In events with several combinatorial 
solutions, only the one with the minimum in the fit is accepted. The main improvement 
relative to the method described in Ref. 151 is the increase in the number of correct assignments 
of b jets to b quarks. This is achieved by applying the kinematic fit twice, sequentially, in each 
event. In the first fit, the top quark mass is fixed to a value of 172.5 GeV. The jet combina¬ 
tion that provides the minimum x^ in the fit is then used as input to the second kinematic fit, 
in which the top quark mass is not fixed, and the solution to this fit is retained. A further 
improvement in the method is to require the ^^'probability of the second kinematic fit to be 
>2%. This criterion is chosen to optimize the fraction of correctly reconstructed signal events, 
without increasing significantly the statistical uncertainty in the data. The efficiency of this 
requirement is about 87% for signal events with the correct jet assignment. As a result, the 
number of correctly reconstructed events is increased by almost a factor of two relative to the 
method used in Ref. IQ, and effects from migration of events across bins, which are relevant 
for the measurements of the cross section, are reduced. It has been checked that any possible 
bias in the results that could be introduced by fixing the top quark mass to a specific value in 
the first kinematic fit is within the assigned systematic uncertainty on the dependence of the 
measurement on the top quark mass (cf. Section [5^ . 

The dilepton channels use an algebraic kinematic reconstruction method |j5lll5|. The only un¬ 
knowns are the three-momenta of the two neutrinos, which are reconstructed imposing the fol¬ 
lowing kinematic constraints: px conservation in the event; the W bosons, and top quark and 
antiquark masses. In contrast to the method of Ref. |5|, the top quark mass is fixed to a value of 
172.5 GeV. Each suitable pair of b jet candidates in the event, and both possible assignments of 
these two jets to the two selected leptons, are considered in the kinematic reconstruction. Gom- 
binations with two b-tagged jets are preferred to using single b-tagged jets. In the new method, 
events are reconstructed 100 times, each time randomly smearing the measured energies and 
directions of the reconstructed lepton and b jet candidates by their respective detector resolu¬ 
tions. This smearing recovers events that yielded no solution of the equations for the neutrino 
momenta, because of measurement fluctuations. The equations for the neutrino momenta can 
have up to four solutions. For a given smearing, the solution is identified by the one yielding 
the smallest invariant mass of the tt system. For each solution, a weight is calculated based 
on the expected true lepton-b-jet invariant mass spectrum. The weights are summed over the 
100 reconstruction attempts, and the kinematic quantities associated to the top quark and an¬ 
tiquark are calculated as a weighted average. Finally, the two jet and lepton-jet assignments 
that yield the maximum sum of weights are chosen for analysis. It has been checked that any 
bias introduced through the use of the lepton-b-jet and tt invariant masses is negligible. This 
method yields on average a reconstruction efficiency of ^9A7o, which is 6% higher than the one 
described in Ref. 0, and reduces systematic migration effects. 
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5 Systematic uncertainties 


Distributions of the top quark or antiquark and tt kinematic observables (the transverse mo¬ 
menta pj, and the rapidities i/t and j/jj) are presented in Figs. and for the ^+jets and 
dilepton channels, respectively. The hatched regions correspond to the shape uncertainties for 
the signal and background (cf. Section|^, and are dominated by the former. The lower panel in 
each plot also shows the ratio of data relative to the simulated signal and background samples. 

In general, the data are reasonably well described by the simulation within the uncertainties. 
For both channels, the measured pj distributions, in particular pj, are somewhat softer than 
the simulated distributions: the data lie above the simulation for pj < 60(65) GeV in the £+jets 
(dilepton) channels, while they lie below for pj > 200 GeV. This pattern was also observed 
at 7TeV IHI. To ensure that the results presented in Section]^ are not affected by such small 
remaining differences between data and simulation, the analysis has been repeated in different 
kinematic regions, with different selection requirements, and after scaling the top quark pj 
spectrum in simulation to match the data. However, the impact on the measurement of the 
cross sections is negligible. 

Following the event selection described in Section |4^ and the kinematic reconstruction of the 
tf system, the main contributions to the background in the £+jets channels arise from tf de¬ 
cays into channel other than f'+jets (including tf decays into t leptons originating from the 
primary interaction) and single top quark events. The contribution from W+jets and QGD 
multijet events are well suppressed after the b-tagging requirement, while other tf events are 
somewhat reduced after the ;y^-probability requirement. A total of 24 927 events are found in 
the e+jets channel and 26 843 events in the p+jets channel. The contribution from tf signal to 
the final event sample is 89.0%. The remaining fraction of events contains 7.3% tf decays other 
than the f+jets channels, 2.4% single top quark events, 0.9% W+jets and tf +Z/W /7 events, and 
negligible fractions of Z+jets, diboson, and QGD multijet events. All background contributions 
are determined from simulation. 

In the dilepton channels, 10 678 events are found in the e^e^ channel, 14403 in the chan¬ 
nel, and 39 640 in the e^p^ channel. Only tf events containing at least two leptons (electrons or 
muons) from W decays in the final state are considered as signal, and constitute 79.0% of the fi¬ 
nal event sample. All other tf candidate events, specifically those originating from decays via r 
leptons, are considered as background and amount to 13.3% of the final event sample. The frac¬ 
tion of Z+jets events is found to be 2.4%. This background, which is dominant to the e+e^ and 
p+p^ channels, is estimated from data using the number of events observed within the Z-peak 
region (which is removed from the candidate sample), and a correction needed for non-Z+jets 
backgrounds in this same control region is obtained from data in the e^p^ channel Il5l l4^ . 
Other sources of background, including single top quark production (3.4%), tf +Z/W /7 pro¬ 
duction ( 1 %), the contribution arising from misidentified or genuine leptons within jets ( 0 . 6 %), 
or diboson events (0.3%), are estimated from simulation. 


5 Systematic uncertainties 

The measurement is affected by systematic uncertainties that originate from detector effects 
and from theoretical assumptions. Each source of systematic uncertainty is assessed individ¬ 
ually by changing the corresponding efficiency, resolution, or scale by its uncertainty, using 
a prescription similar to the one followed in Ref. ||5l. For each change made, the measured 
normalized differential cross section is recalculated, and the difference of the changed result 
relative to its nominal value in each bin is taken as the systematic uncertainty. The overall un¬ 
certainty on the measurement is obtained by adding all the contributions in quadrature, and is 
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the kinematic reconstruction in the £+jets charmels. The top row shows the pj, and the bot¬ 
tom row shows the rapidities. The QCD multijet background is negligible and not shown. 
The hatched regions correspond to the shape uncertainties for the signal and backgrounds (cf. 
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5 Systematic uncertainties 
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Figure 4: Distribution of top quark or antiquark (left) and tt (right) quantities as obtained from 
the kinematic reconstruction in the dilepton channels. The top row shows the pj, and the 
bottom row shows the rapidities. The QCD multijet background is negligible and not shown. 
The Z/ 7 *+jets background is determined from data IHlIlll. The hatched regions correspond to 
the shape uncertainties for the signal and backgrounds (cf. Section]^. The lower part of each 
plot shows the ratio of data to the predictions. 
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of the order of 3-10%, depending on the observable and the bin. A detailed description of this 
is given in Sections [5d] and 5.2 The typical representative values of the systematic uncertainties 
in the normalized differential cross sections are summarized in Table [T] 


5.1 Experimental uncertainties 


The efficiencies of the single-electron and single-muon triggers in the ^+jets channels are de¬ 
termined using the "tag-and-probe" method of Ref. using Z boson event samples. Scale 
factors close to unity within a few percent are extracted to account for the observed dependence 
on the t] and px of the lepton. The lepton identification and isolation efficiencies for the £+jets 
channels obtained with the tag-and-probe method agree well between data and simulation, so 
that the applied corrections are very close to unity The systematic uncertainties are determined 
by shape-dependent changes in trigger and selection efficiencies by their uncertainties. Lepton 
trigger efficiencies in the dilepton channels are measured using triggers that are only weakly 
correlated to the dilepton triggers used in the analysis. A dependence on p of a few percent is 
observed, and scale factors are extracted. The lepton identification and isolation uncertainties 
in the dilepton channels are also determined using the tag-and-probe method, and are again 
found to be described very well by the simulation for both electrons and muons. The overall 
difference between data and simulation in bins of t] and px is estimated to be < 2 % for electrons, 
and scale factors for muons are found to be close to unity within 1 . 0 %. 

The uncertainty due to the limited knowledge of the jet energy scale is determined by changes 
implemented in jet energy in bins of px and rj [I38II . The uncertainty due to the limited accu¬ 
racy of the jet energy resolution (JER) is determined by changing the simulated JER by ±ln in 
different rj regions II38II . 

The uncertainty in b-tagging efficiency is determined by taking the maximum change in the 
shape of px and t] b jet distributions obtained by changing the scale factors. This is achieved 
by dividing the b jet distributions in px and f] into two bins at the median of the respective 
distributions. These correspond to px = 65GeV, and |p| = 0.7 and 0.75 for the £+jets and 
dilepton channels, respectively The b-tagging scale factors for b jets in the first bin are scaled 
up by half of the uncertainties quoted in Ref. Il40l . while those in the second bin are scaled 
down, and vice versa, so that a maximum variation is assumed and the difference between the 
scale factors in the two bins reflects the full uncertainty The changes are made separately in 
the Px and rj distributions, and independently for heavy-flavour (b and c) and light (s, u, d, and 
gluon) jets, assuming that they are all uncorrelated. 


The uncertainty in background normalization is determined by changing the background yields. 
In the £+jets channels, the background normalization for the diboson, QCD multijet, W+jets, 
and Z+jets samples is conservatively varied by ±50% O, since these backgrounds, being very 
small, are determined from simulation rather than from data. The normalization of the tf 
+Z/W /7 samples is changed by ±30%. Eor the single top quark sample, the uncertainty is 
covered by changing the normalization by ±30%, and the kinematic scales of the event process 
(renormalization and factorization scales) as described in Section 5.2 In the e^e^ and p^p~ 
channels, the dominant background from Z+jets determined from data |j5l |1^ is changed in 
normalization by ±30%. In addition, changes in the background contributions from single top 
quark, diboson, QCD multijet, tt +Z/W/ 7 , and W+jets events of ±30% are used in dilepton 
channels 151. 


The kinematic reconstruction of top quarks is well described by the simulation, and the result¬ 
ing uncertainties are small. In the case of the £+jets analysis, the uncertainty of the kinematic 
fit is included in the changes in jet energy scales and resolutions, and in the uncertainty on the 
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6 Normalized differential cross sections 


dependence on the top quark mass (cf. Section [5^ . In the dilepton analysis, the bin-to-bin un¬ 
certainty is determined from the small remaining difference in efficiency between simulation 
and data. 

The pileup model estimates the mean number of additional pp interactions to be about 20 
events per bunch crossing for the analyzed data. This estimate is based on the total inelastic 
proton-proton cross section, which is determined to be 69.4 mb following the measurement 
described in Ref. HTl . The systematic uncertainty is determined by changing this cross section 
within its uncertainty of ±5%. 


5.2 Uncertainties in modeiiing 


The impact of theoretical assumptions on the measurement is determined, as indicated pre¬ 
viously, by repeating the analysis and replacing the standard MadGraph tt simulation by 
dedicated simulation samples with altered parameters. 


The uncertainty in modelling of the hard-production process is assessed through changes in 
the renormalization and factorization scales in the MadGraph sample by factors of two and 
0.5 relative to their common nominal value, which is set to the Q of the hard process. In Mad¬ 
Graph, Q is defined by = m\ + "^Pj, where the sum is over all additional final state partons 
in the matrix element. The impact of the choice of the scale that separates the description of jet 
production through matrix elements (ME) or parton shower (PS) in MadGraph is studied by 
changing its reference value of 20 GeV to 40 GeV and to 10 GeV. In the £+jets channels, changes 
in the renormalization and factorization scales are also applied to single top quark events to 
determine an uncertainty on the shape of this background contribution. The dependence of 
the measurement on the top quark mass is also estimated from dedicated MadGraph simula¬ 
tion samples in which the top quark mass is changed by ±1 GeV relative to the value used in 
the default simulation. The uncertainty from hadronization and parton showering is assessed 
by comparing the results obtained from samples simulated with POWHEG and MC@NLO in¬ 
terfaced with PYTHIA6 and HERWIG6, respectively. The uncertainty from the choice of PDF is 
determined by reweighting the sample of simulated tf signal events according to the 52 GTIO 
PDF error sets Il24l , at a 90% confidence level. The maximum variation is taken as uncertainty 


As mentioned in Sections 4.2 and |4.3} the effect of scaling the top quark pj spectrum in sim¬ 
ulation to match the data has negligible impact on the measured cross sections, therefore no 
systematic uncertainty is taken into account for this effect. 


6 Normalized differential cross sections 

The normalized tf cross section in each bin i of each observable X is determined as a function 
of the kinematic properties of the leptons, the lepton pair, the b jets, the b jet system, the top 
quarks, and the tf system through the relation Uj: 

1 daj _ 1 Xj 

adX EiXiAf 

where xi represents the number of signal events measured in data in bin i after background 
subtraction and corrected for detector efficiencies, acceptances, and migrations, and Af is the 
bin width. The differential cross section is normalized by the sum of x, over all bins, as indi¬ 
cated in Eq. Q. The integrated luminosity is omitted, as it cancels in the ratio. Because of the 
normalization, sources of systematic uncertainty that are correlated across all bins of the mea¬ 
surement, e.g. the uncertainty in the integrated luminosity, also cancel. The contribution to the 
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Table 1: Breakdown of typical systematic uncertainties for the normalized differential cross sec¬ 
tions. The uncertainty in the jet-parton matching threshold is indicated as "ME-PS threshold"; 
"PS" refers to "parton shower". The medians of the distribution of uncertainties over all bins of 
the measurement are quoted. For the £+jets channels, the background from Z+jets is negligible 
and included in the "Background (all other)" category. 


Relative systematic uncertainty (%) 


Source 

Fepton 

and b jet observables 

Top quark and tt observables 


£+jets 

dileptons 

£+jets 

dileptons 

Trigger eff. & lepton selec. 

0.1 

0.1 

0.1 

0.1 

Jet energy scale 

2.3 

0.4 

1.6 

0.8 

Jet energy resolution 

0.4 

0.2 

0.5 

0.3 

Background (Z+jets) 

— 

0.2 

— 

0.1 

Background (all other) 

0.9 

0.4 

0.7 

0.4 

b tagging 

0.7 

0.1 

0.6 

0.2 

Kinematic reconstruction 

— 

<0.1 

— 

<0.1 

Pileup 

0.2 

0.1 

0.3 

0.1 

Fact./renorm, scale 

1.1 

0.7 

1.8 

1.2 

ME-PS threshold 

0.8 

0.5 

1.3 

0.8 

Hadronization & PS 

2.7 

1.4 

1.9 

1.1 

Top quark mass 

1.5 

0.6 

1.0 

0.7 

PDF choice 

0.1 

0.2 

0.1 

0.5 


background from other tt decays is taken into account, after subtracting all other background 
components, by correcting the number of signal events in data using the expected signal frac¬ 
tion. The expected signal fraction is defined as the ratio of the number of selected tf signal 
events to the total number of selected tt events (i.e. signal and all other tt events) in simula¬ 
tion. This procedure avoids the dependence on the total inclusive tf cross section used in the 
normali z ation of the simulated signal sample. 

Effects from trigger and detector efficiencies and resolutions leading to the migration of events 
across bin boundaries, and therefore to statistical correlations among neighbouring bins, are 
corrected by using a regularized unfolding method ||5lll8lll3. For each measured distribution, 
a response matrix is defined that accounts for migrations and efficiencies using the simulated 
MadGraph+pythia6 tf signal sample. The generalized inverse of the response matrix is used 
to obtain the unfolded distribution from the measured distribution by applying a mini¬ 
mization technique. A smoothing prescription (regularization) is applied to prevent large un¬ 
physical fluctuations that can be introduced when directly inverting the response matrix. The 
strength of the regularization is determined and optimized individually for each distribution 
using the averaged global correlation method [50l. To keep the bin-to-bin migrations small, the 
widths of bins in the measurement are chosen according to their purity (ratio of the number 
of events generated and reconstructed in a particular bin to the total number of events recon¬ 
structed in that bin; this quantity is sensitive to migrations into the bin) and stability (ratio of 
the number of events generated and reconstructed in a particular bin to the number of events 
generated in that bin; this is sensitive to migrations out of the bin). The purity and stability 
of the bins in this analysis are typically 60% or larger, mainly due to the improvements in the 
kinematic reconstruction methods discussed in Section 1431 

The performance of the unfolding procedure is tested for possible biases from the choice of 
the input model (the MadGraph+pythia6 tf signal simulation). It is verified that, either by 
reweighting the signal simulation or injecting a resonant tf peak into the simulation of the sig- 
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6 Normalized differential cross sections 


nal, the unfolding procedure based on the nominal response matrices still recovers these altered 
shapes within statistical uncertainties. Moreover, tt samples simulated with powheg+pythia6 
and mc@NLO+herwig6 are used to obtain the response matrices applied in the unfolding 
when determining the systematic uncertainties of the model (cf. Section [5^ . Therefore, pos¬ 
sible effects from the unfolding procedure are already taken into account in the systematic 
uncertainties. The unfolded results are found to be consistent with those obtained using other 
regularization techniques |49|. 

The measurement of the normalized differential cross sections proceeds as follows. For each 
kinematic distribution, the event yields in the separate channels are added together, the back¬ 
ground is subtracted, and the unfolding is performed. It is verified that the measurements in 
separate channels yield results consistent within their uncertainties. The systematic uncertain¬ 
ties in each bin are determined from the changes in the combined cross sections. This requires 
the full analysis to be repeated for every systematic change, and the difference relative to the 
nominal combined value is taken as the systematic uncertainty for each bin of each observ¬ 
able. This method therefore takes into account the correlation among systematic uncertainties 
in different charmels and bins. 

The normalized differential cross sections of leptons and b jets are unfolded to the particle 
level and determined in a fiducial phase space defined by the kinematic and geometric re¬ 
gion in which the final-state leptons and jets are produced within the detector acceptance 
(cf. Section [6d] |. This minimizes model uncertainties from the extrapolation of the measure¬ 
ment outside of the experimentally well-described regions of phase space. In addition, the 
top quark and tt-system quantities are unfolded to the parton level and presented in the full 
phase space (cf. Section [6!^ to provide easier comparisons with recent QCD calculations. The 
measurements are compared to predictions from MadGraph+PYTHIA6, P0WHEG+PYTHIA6, 
POWHEG+HERWIG6, and MC@NL0+HERWIG6. The top quark and tf results are also compared 
to the latest calculations at NLO+NNLL [Til [HI and approximate NNLO IT^ precision, when 
available. 

In addition to the measurements discussed in Ref. [Sj, results for the pj and invariant mass of 
the b jet pair, the pj of the top quarks or antiquarks in the tf rest frame, the pj of the highest 
(leading) and second-highest (trailing) pj of the top quark or antiquark, and the difference in 
the azimuthal angle between the top quark and antiquark are also presented. 

All values of normalized differential cross sections, including bin boundaries, are provided in 
tables in Appendix]^ 


6.1 Lepton and b jet differential cross sections 


The normalized differential tf cross section as a function of the lepton and b jet kinematic prop¬ 
erties is measured at the particle level, where the objects are defined as follows. Leptons from W 
boson decays are defined after final-state radiation. A jet is defined at the particle level, follow¬ 
ing a procedure similar to that described in Section 4.1 for reconstructed jets, by applying the 
anti-kx clustering algorithm with a distance parameter of 0.5 to all stable particles (excluding 
the decay products from W boson decays into ei/, pv, and final states with leptonic t decays). 
A jet is defined as a b jet if it contains any of the decay products of a B hadron. Only the two 
b jets of highest px originating from different B hadrons are considered as arising from the top 
quark decays. 


The measurements are presented in a fiducial phase space defined by geometric and kinematic 
requirements on these particle-level objects as follows. The charged leptons from the W boson 
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decays must have |?y| < 2.1 and pj > 33GeV in the £+jets channels, and |p| < 2.4 and pj > 
20 GeV in the dilepton channels. Exactly one and two leptons are required, respectively, in 
the £+jets and the dilepton channels. At least four jets with |p| < 2.4 and px > 30 GeV, two 
of which are b jets, are required in the £+jets channels. In the dilepton channels, both b jets 
from the top quark decays must satisfy |p| <2.4 and px > 30 GeV. The fiducial particle-level 
corrections are determined using simulated tt events that fulfill these requirements; all other tf 
events are classified as background and are removed. 

Figure presents the normalized differential cross section in the ^+jets channels as a func¬ 
tion of the lepton transverse momentum pj and pseudorapidity rjg. The distributions of the 
transverse momentum of the b jets pj and their pseudorapidity pb ^re given in Fig. to¬ 
gether with the transverse momentum p^ and invariant mass of the b jet pair. Also 
shown are predictions from MadGraph+PYTHIA6, POWHEG+PYTHIA6, POWHEG+HERWIG6, 
and mc@NLO+herwig6. The lower panel in each plot shows the ratio of each of these predic¬ 
tions to data, in order to quantify their level of agreement relative to data. 


CMS, 19.7fb'at fs = 8TeV 


CMS, fs = 8TeV 



[GeV] 



Figure 5: Normalized differential tt production cross section in the £+jets channels as a function 
of the pj (left) and (right) of the charged lepton. The superscript '£' refers to both £+ and 
. The data points are placed at the midpoint of the bins. The inner (outer) error bars indicate 
the statistical (combined statistical and systematic) uncertainties. The measurements are com¬ 
pared to predictions from MadGraph+PYTHIA6, POWHEG+PYTHIA6, POWHEG+HERWIG6, 
and mc@NLO+herwig6. The lower part of each plot shows the ratio of the predictions to 
data. 


Figure [^presents the normalized differential cross sections for the dilepton charmels: the trans¬ 
verse momentum pj and the pseudorapidity p^ of the leptons, and the transverse momentum 
p^^ and the invariant mass ni£+£- of the lepton pair. The distributions in the transverse mo¬ 
mentum of the b jets pj and their pseudorapidity pb are shown in Fig. together with the 
transverse momentum p^ and invariant mass of the b jet pair. Predictions from Mad- 
Graph+PYTHIA6, POWHEG+PYTHIA6, POWHEG+HERWIG6, and MC@NL0+HERWIG6 are also 
presented for comparison. 

In general, none of the examined predictions provides an accurate description of data for all 
measured lepton and b jet distributions. A steeper px spectrum is observed in data for the 
lepton and the b jet distributions compared to the predictions in both decay channels, which 
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6 Normalized differential cross sections 
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Figure 6: Normalized differential tf production cross section in the £+jets channels as a function 
of the pj (top left) and r]\, (top right) of the b jets, and the p^ (bottom left) and (bottom right) 

of the b jet pair. The superscript 'b' refers to both b and b jets. The data points are placed at the 
midpoint of the bins. The inner (outer) error bars indicate the statistical (combined statistical 
and systematic) uncertainties. The measurements are compared to predictions from Mad- 
GRAPH+PYTHIA6, POWHEG+PYTHIA6, POWHEG+HERWIG6, and MC@NL0+HERWIG6. The 
lower part of each plot shows the ratio of the predictions to data. 
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is best described by powheg+herwig6. The lepton pj in data is above the predictions for 
Pj < 40 GeV, while it is below for pj > 100 GeV. A similar behaviour is observed for p^^ , 
Pj, and p^. The ni^+i- distribution in data is below all predictions for m£+£- > 30 GeV. Worse 
agreement is found for POWHEG+PYTHIA6. The p distributions in data are described by the 
predictions within the experimental uncertainties. The pb distributions are slightly less cen¬ 
tral in data than in the predictions, and are worse described by MadGraph+pythia6. The 
remaining distributions are described by the predictions within experimental uncertainties. 

6.2 Top quark and tt differential cross sections 

The normalized differential tf cross section as a function of the kinematic properties of the top 
quarks and the tt system is defined with respect to the top quarks or antiquarks before the 
decay (parton level) and after QGD radiation, and extrapolated to the full phase space using 
the MadGraph+pythia 6 prediction for the £+jets and dilepton charmels. 

In Figs.|^to[^ the following distributions are presented for the ^+jets channels: the transverse 
momentum and the rapidity i/t of the top quarks or antiquarks, the transverse momen¬ 
tum pj of the top quarks or antiquarks in the tf rest frame, the difference in the azimuthal 
angle between the top quark and antiquark A(p(t,f), the transverse momentum of the leading 
(Pj) and trailing (p^) top quark or anti quark, and the transverse momentum pj, the rapidity 
Ptf, and the invariant mass m^i of the tf system. The data are compared to predictions from 
MADGRAPH+PYTHIA6, POWHEG+PYTHIA6, POWHEG+HERWIG6, and MC@NL0+HERWIG6. In 
addition, the approximate NNLO calculation fl^ is also shown for the top quark pj and rapid¬ 
ity results, while the m££ and the pj distributions are compared to the NLO+NNLL predictions 
from Refs. I14| and HTSl , respectively. Figures p^] - [l4| show the corresponding distributions in 
the dilepton channels. The lower panel in each plot also shows the ratio of each prediction 
relative to data. 

In general, the powheg+herwig6 prediction provides a good description of data for all mea¬ 
sured distributions. The shape of the top quark pj spectrum is softer in data than in the predic¬ 
tions from MadGraph+PYTHIA6, POWHEG+PYTHIA6, and MC@NL0+HERWIG6 in both chan¬ 
nels. The data lie above the predictions for pj < 60 (65) GeV in the f+jets (dilepton) channels, 
while they lie below for pj > 200 GeV. This effect was also observed at 7TeV |5J. The disagree¬ 
ment between data and predictions in the tail of the distributions is also observed in a mea¬ 
surement by the ATLAS Gollaboration ||6|. In contrast, the prediction from P0WHEG+HERWIG6 
and the approximate NNLO calculation provide a better description of the data, as they pre¬ 
dict a slightly softer top quark pj distribution than the three other simulations. The difference 
between the P0WHEG+PYTHIA6 and P0WHEG+HERWIG6 distributions is attributed to differ¬ 
ent treatment of the hardest initial state radiation in PYTHIA6 and HERWIG6. The same pat¬ 
tern is observed for pj, indicating that the softer spectrum in data is not caused by the boost 
of the tt system. It is also present in the pj, and particularly, in the pj distributions. For 
all these distributions, the powheg+herwig6 prediction provides a better description of the 
data. The difference in the shape of the top quark pj spectrum between data and simulation 
is observed consistently in the analyses using different event selection requirements or differ¬ 
ent pileup conditions. The yt distribution is found to be slightly less central in data than in 
the predictions, particularly in the case of MadGraph+pythia6 and the approximate NNLO 
calculation, which are more central than the other predictions. On the contrary, is more cen¬ 
tral in data, and it is slightly better described by MadGraph+pythia6. The distribution 
in data tends to be lower than the predictions for large values, and is better described by 
MadGraph+pythia6 and powheg+herwig6. The p^ spectrum is well described by all the 
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CMS, (s= SleV CMS, 19.7fb 'at 8TeV 




Figure 7: Normalized differential ft production cross section in the dilepton channels as a func¬ 
tion of the pj (top left) and rji (top right) of the charged leptons, and the p^^ (bottom left) 
and m£+£- (bottom right) of the lepton pair. The superscript refers to both and The 
data points are placed at the midpoint of the bins. The irmer (outer) error bars indicate the 
statistical (combined statistical and systematic) uncertainties. The measurements are com¬ 
pared to predictions from MadGraph+PYTHIA6, POWHEG+PYTHIA6, POWHEG+HERWIG6, 
and mc@NLO+herwig6. The lower part of each plot shows the ratio of the predictions to 
data. 
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Figure 8: Normalized differential tt production cross section in the dilepton channels as a 
function of the p\ (top left) and (top right) of the b jets, and the (bottom left) and 
?«bb (bottom right) of the b jet pair. The superscript 'b' refers to both b and b jets. The 
data points are placed at the midpoint of the bins. The mner (outer) error bars indicate the 
statistical (combined statistical and systematic) uncertainties. The measurements are com¬ 
pared to predictions from MadGraph-i-PYTHIA6, POWHEG-i-PYTHIA6, POWHEG-i-HERWIG6, 
and mc@NLO+herwig6. The lower part of each plot shows the ratio of the predictions to 
data. 
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7 Summary 


considered predictions, except for the NLO+NNLL calculation, which fails to describe the data 
for all pj values. 


The results from the £+jets and dilepton channels are compared to each other in Figs. 15 to 17 


This is only feasible for the top quark and tt quantities, since they are measured in the same 
phase space (i.e. the full parton level phase space) for both channels. The results are presented 
relative to the MadGraph+pythia6 prediction to highlight the level of agreement between 
data and the default tt simulation. To facilitate the comparison of measurements that are per¬ 
formed using different size and number of bins, a horizontal bin-centre correction is applied 
to all data points from both channels. In each bin, the measured data points are presented at 
the horizontal position in the bin where the predicted bin-averaged cross section equals the 
cross section of the unbinned MadGraph+PYTHIA6 calculation (cf. Ell), which is common 
for both channels. The data are also compared to the predictions from powheg+pythia6, 
POWHEG+HERWIG6, MC@NL0+HERWIG6 relative to MadGraph+PYTHIA6. The results are 
consistent between the channels for all quantities, in particular, for all measurements related to 
the top quark pj distribution. The softer spectrum in data relative to MadGraph+pythia6 is 
clearly visible. 


In addition, a comparison between results obtained at \/s = 7 |51 and 8 TeV is also performed 
for both the ^+jets and dilepton channels, and presented in Figs. 18 and 19 for pj, yt, pj, i/^. 


and JMjf. Since the fiducial phase space definition for the normalized differential cross sections 
is also different for each value of \/s, the comparison is again only possible for top quark and 
tt quantities. The measurements are presented relative to the corresponding default Mad- 
Graph+PYTHIA6 predictions at 7 and 8 TeV. A horizontal bin-centre correction with respect to 
the MadGraph+pythia6 predictions is applied to all data points from both channels and 
values. The results are consistent between the channels for all quantities, both at 7 and 8 TeV. 
The uncertainties in almost all bins of the distributions are reduced for the 8 TeV results relative 
to 7 TeV, mainly due to the improvements discussed in Section [43} The softer pj in data relative 
to MadGraph+pythia6 is also visible at 7 TeV. 


7 Summary 

First measurements are presented of normalized differential tf production cross sections in 
pp collisions at ^/s = 8 TeV. The measurements are performed with the GMS detector in the 
£+jets {i = e or p) and dilepton (e+e^, and e^p^) tt decay channels. The normalized 

tf cross section is measured as a function of the transverse momentum, rapidity, and invariant 
mass of the final-state leptons and b jets in the fiducial phase space, and the top quarks and tt 
system in the full phase space. The measurements in the different decay channels are in agree¬ 
ment with each other. In general, the data are in agreement with standard model predictions 
up to approximate NNLO precision. Among the examined predictions, P0WHEG+HERWIG6 
provides the best overall description of the data. However, the pj spectrum in data for 
leptons, jets, and top quarks is softer than expected, particularly for MadGraph+pythia6, 
POWHEG+PYTHIA6, and MC@NL0+HERWIG6. The calculation at approximate NNLO precision 
also provides a good description of the top quark pj spectrum. The distribution in data 
tends to be lower than the predictions for large m^l values. The pj spectrum is well described 
by all the considered predictions, except for the NLO+NNLL calculation, which fails to de¬ 
scribe the data for all pj values. The results show the same behaviour as the corresponding 
GMS measurements at ^/s = 7 TeV. 
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Figure 9: Normalized differential tt production cross section in the £+jets channels as a function 
of the pj (top left), the tf rest frame pj (top right), and the rapidity x/t (bottom left) of the 
top quarks or antiquarks, and the difference in the azimuthal angle between the top quark 
and the antiquark A^(t,f) (bottom right). The data points are placed at the midpoint of the 
bins. The inner (outer) error bars indicate the statistical (combined statistical and systematic) 
uncertainties. The measurements are compared to predictions from MadGraph+pythia6, 
POWHEG+PYTHIA6, POWHEG+HERWIG6, MC@NL0+HERWIG6, and to approximate NNLO HU 
calculations, when available. The lower part of each plot shows the ratio of the predictions to 
data. 
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7 Summary 



Figure 10: Normalized differential tf production cross section in the £+jets channels as a 
function of the pj of the leading (left) and trailing (right) top quarks or antiquarks. The 
data points are placed at the midpoint of the bins. The irmer (outer) error bars indicate the 
statistical (combined statistical and systematic) uncertainties. The measurements are com¬ 
pared to predictions from MadGraph+PYTHIA6, POWHEG+PYTHIA6, POWHEG+HERWIG6, 
and mc@NLO+herwig6. The lower part of each plot shows the ratio of the predictions to 
data. 
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Figure 11: Normalized differential ft production cross section in the £+jets channels as a func¬ 
tion of the pj (top left), (top right), and (bottom) of the tt system. The data points are 
placed at the midpoint of the bins. The inner (outer) error bars indicate the statistical (combined 
statistical and systematic) uncertainties. The measurements are compared to predictions from 
MADGRAPH+PYTHIA6, POWHEG+PYTHIA6, POWHEG+HERWIG6, MC@NL0+HERWIG6, and to 
NLO+NNLL 1(141 [l5l calculations, when available. The lower part of each plot shows the ratio 
of the predictions to data. 
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Figure 12: Normalized differential tf production cross section in the dilepton channels as a 
function of the pj (top left), the tf rest frame pj (top right), and the rapidity i/t (bottom left) of 
the top quarks or antiquarks, and the difference in the azimuthal angle between the top quark 
and the antiquark A^(t,f) (bottom right). The data points are placed at the midpoint of the 
bins. The inner (outer) error bars indicate the statistical (combined statistical and systematic) 
uncertainties. The measurements are compared to predictions from MadGraph+pythia6, 
POWHEG+PYTHIA6, POWHEG+HERWIG6, MC@NL0+HERWIG6, and to approximate NNLO HU 
calculations, when available. The lower part of each plot shows the ratio of the predictions to 
data. 
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Figure 13: Normalized differential ft production cross section in the dilepton channels as a 
function of the pj of the leading (left) and trailing (right) top quarks or antiquarks. The 
data points are placed at the midpoint of the bins. The irmer (outer) error bars indicate the 
statistical (combined statistical and systematic) uncertainties. The measurements are com¬ 
pared to predictions from MadGraph+PYTHIA6, POWHEG+PYTHIA6, POWHEG+HERWIG6, 
and mc@NLO+herwig6. The lower part of each plot shows the ratio of the predictions to 
data. 
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Figure 14: Normalized differential ft production cross section in the dilepton channels as a 
function of the pj (top left), (top right), and (bottom) of the tt system. The data points are 
placed at the midpoint of the bins. The inner (outer) error bars indicate the statistical (combined 
statistical and systematic) uncertainties. The measurements are compared to predictions from 
MADGRAPH-I-PYTHIA6, POWHEG-I-PYTHIA6, POWHEG-I-HERWIG6, MC@NL0-I-HERWIG6, and to 
NLO-lNNLL IfT^llSl calculations, when available. The lower part of each plot shows the ratio 
of the predictions to data. 
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Figure 15: Comparison of normalized differential tf production cross section in the dilepton 
and £+jets channels as a function of the pj (top left), the tf rest frame pj (top right), and 
the rapidity x/t (bottom left) of the top quarks or antiquarks, and the difference in the az¬ 
imuthal angle between the top quark and the antiquark A^(t,f) (bottom right). The measure¬ 
ments are presented relative to the MadGraph+pythia6 prediction. A horizontal bin-centre 
correction is applied to all data points (cf. Section |6.2| |. The inner (outer) error bars indi¬ 
cate the statistical (combined statistical and systematic) uncertainties. The predictions from 
POWHEG+PYTHIA6, POWHEG+HERWIG6, and MC@NL0+HERWIG6, also presented relative to 
MadGraph+PYTHIA6, are shown for comparison. 
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Figure 16: Comparison of normalized differential tf production cross section in the dilepton 
and £+jets charmels as a function of the pj of the leading (left) and trailing (right) top quarks or 
antiquarks. The measurements are presented relative to the MadGraph+pythia6 prediction. 
A horizontal bin-centre correction is applied to all data points (cf. Section |6.2| . The inner 
(outer) error bars indicate the statistical (combined statistical and systematic) uncertainties. 
The predictions from POWHEG+PYTHIA6, POWHEG+HERWIG6, and MC@NL0+HERWIG6, also 
presented relative to MadGraph+pythia6, are shown for comparison. 
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Figure 17: Comparison of normalized differential ft production cross section in the dilepton 
and ^+jets channels as a function of the pj (top left), (top right), and (bottom) of the 
tf system. The measurements are presented relative to the MadGraph+pythia6 prediction. 
A horizontal bin-centre correction is applied to all data points (cf. Section |6.2| . The inner 
(outer) error bars indicate the statistical (combined statistical and systematic) uncertainties. 
The predictions from POWHEG+PYTHIA6, POWHEG+HERWIG6, and MC@NL0+HERWIG6, also 
presented relative to MadGrapei+pythia6, are shown for comparison. For better visibility, 
data points with identical bin centres (cf. Tables A.6 and A.lOl are shifted horizontally by a 
negligible amount. 


















30 


7 Summary 


to 

cc 


1.5r 

1 . 4 : 


ehcTi.Sr 

T3| -o 

1.2 ri 


1.1: 


0.9 r 
0.8 r 

0.7 r 

0.6 r 

0.5, 


CMS, 5.0/19.7 fb'' at \fs = 7/8 TeV 

" ' I " " I " " I " " I " " I' " ' I' " ' I' " ' I'' 

Data or theory / MadGraph+Pythiah 
• e/p + Jets (8 TeV) 

A e/p + Jets (7 TeV) 

■ Dilepton (8 TeV) 

T Dilepton (7 TeV) 


^ ' 


0 50 100 150 200 250 300 350 400 450 500 


p; [GeV] 


CMS, 5.0/19.7 fb-' at ifs = 7/8 TeV 



Figure 18: Comparison of normalized differential ft production cross section in the dilepton 
and £+jets channels at 7 TeV [jSlI and 8 TeV, as a function of the pj (left) and rapidity yt (right) of 
the top quarks or antiquarks. The measurements are presented relative to the corresponding 
MadGraph+pythia6 predictions. A horizontal bin-centre correction is applied to all data 
points (cf. Section [6^ . The inner (outer) error bars indicate the statistical (combined statistical 
and systematic) uncertainties. For better visibility, data points with identical bin centres (cf. 
and A.10[| are shifted horizontally by a negligible amount. 


Tables 


A.6 
















31 


CMS, 5.0/19.7 fb-’ at (s = 7/8 TeV CMS, 5.0/19.7 fb’' at fs = 7/8 TeV 



CMS, 5.0/19.7 fb'^ attfs = 7/8 TeV 



Figure 19: Comparison of normalized differential ft production cross section in the dilepton 
and £+jets charmels at 7 TeV |0 and 8 TeV, as a function of the pj (top left), (top right), and 
(bottom) of the tf system. The measurements are presented relative to the corresponding 
MadGraph+pythia6 predictions. A horizontal bin-centre correction is applied to all data 
points (cf. Section [6^ . The inner (outer) error bars indicate the statistical (combined statistical 
and systematic) uncertainties. For better visibility, data points with identical bin centres (cf. 
Tables A.9 and A.12[| are shifted horizontally by a negligible amount. 
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A Values of the normalized differential cross sections 

In Tables |A]^ to |A]T^ the results for the normalized differential distributions are summarized. 
For each distribution, the result of the measurement, together with the bin range, as well as the 
statistical, systematic, and total uncertainties, are provided. For the top quark or antiquark and 
tt quantities, the bin centres corrected according to the MadGraph+pythia6 prediction (cf. 
Section [6^ are also presented. 

Table A.l: Normalized differential tt cross section in the £+jets channels as a function of the 
charged lepton transverse momentum (pj) and pseudorapidity {rji). The superscript 'i' refers 
to both £+ and The results are presented at particle level in the fiducial phase space. The 
statistical and systematic uncertainties are added in quadrature to yield the total uncertainty 


Pj bin range [GeV] 

l/crdcr/dpj [GeV 

Stat. [%] 

Syst. [%] 

Total [%] 

[30,37) 

1.43 X 10-2 

2.2 

2.4 

3.2 

[37,45) 

2.21 X 10-2 

1.4 

2.8 

3.1 

[45,55) 

1.83 X 10-2 

1.4 

1.7 

2.3 

[55,68) 

1.36 X 10-2 

1.5 

1.4 

2.1 

[68,80) 

9.42 X 10-3 

2.0 

1.5 

2.5 

[80,100) 

5.89 X 10-3 

1.8 

1.2 

2.2 

[100,135) 

2.51 X 10-3 

2.1 

6.3 

6.7 

[135,200) 

5.73 X 10-4 

3.3 

9.5 

10.1 

pi bin range 

l/ada/dpi 

Stat. [%] 

Syst. [%] 

Total [%] 

[-2.1,-1.8) 

9.76 X 10-2 

4.3 

4.3 

6.0 

[-1.8,-1.5) 

1.37 X 10-4 

3.3 

3.1 

4.5 

[-1.5,-1.2) 

1.94 X 10-4 

2.7 

2.4 

3.6 

[-1.2,-0.9) 

2.39 X 10-4 

2.5 

2.8 

3.7 

[-0.9,-0.6) 

3.03 X 10-4 

2.0 

2.3 

3.0 

[-0.6,-0.3) 

3.41 X 10-4 

1.9 

1.4 

2.3 

[-0.3,0.0) 

3.60 X 10-4 

1.9 

1.4 

2.4 

[0.0,0.3) 

3.39 X 10-4 

2.0 

2.0 

2.8 

[0.3,0.6) 

3.35 X 10-4 

1.9 

3.0 

3.6 

[0.6,0.9) 

3.06 X 10-4 

1.9 

4.2 

4.6 

[0.9,1.2) 

2.25 X 10-4 

2.2 

2.2 

3.2 

[1.2,1.5) 

2.01 X 10-4 

2.6 

2.2 

3.4 

[1.5,1.8) 

1.36 X 10-4 

3.3 

2.5 

4.2 

[1.8,2.1) 

9.31 X 10-2 

4.2 

5.6 

7.0 
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Table A.2: Normalized differential ft cross section in the £+jets channels as a function of the b jet 
transverse momentum (pj) and pseudorapidity (pb)/ arid the transverse momentum (p^) and 
the invariant mass of the bb system. The superscript 'h' refers to both b and b jets. The 
results are presented at particle level in the fiducial phase space. The statistical and systematic 
uncertainties are added in quadrature to yield the total uncertainty. 


Pj bin range [GeV] 

l/ada/dpj [GeV ^] 

Stat. [%] 

Syst. [%] 

Total [%] 

[30,48) 

1.40 X 

1.2 

7.0 

7.1 

[48,75) 

1.24 X 

0.9 

1.8 

2.0 

[75,180) 

3.67 X 10-3 

0.7 

4.9 

5.0 

[180,400) 

1.34 X 10-4 

3.3 

12.4 

12.9 

Pb bin range 

l/udcVdpb 

Stat. [%] 

Syst. [%] 

Total [%] 

[-2.4,-1.5) 

1.01 X 10-4 

1.7 

2.9 

3.3 

[-1.5,-1.0) 

2.08 X 10-4 

1.3 

1.5 

2.0 

[-1.0,-0.5) 

2.77 X 10-4 

1.1 

1.2 

1.7 

[-0.5,0.0) 

3.23 X 10-4 

1.1 

2.0 

2.2 

[0.0,0.5) 

3.28 X 10-4 

1.1 

1.9 

2.1 

[0.5,1.0) 

2.89 X 10-4 

1.1 

1.0 

1.5 

[1.0,1.5) 

2.09 X 10-4 

1.4 

2.2 

2.6 

[1.5,2.4) 

1.03 X 10-4 

1.6 

2.8 

3.2 

p^ bin range [GeV] 

l/ada/dp^ [GeV-4] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0,35) 

3.41 X 10-3 

2.4 

4.1 

4.7 

[35,75) 

6.55 X 10-3 

1.7 

3.0 

3.5 

[75,115) 

8.51 X 10-3 

1.5 

3.0 

3.4 

[115,155) 

5.04 X 10-3 

2.0 

6.7 

7.0 

[155,280) 

5.79 X 10-4 

3.8 

11.3 

11.9 

[280,500) 

0.17 X 10-4 

14.1 

23.3 

27.2 

bin range [GeV] 

l/ada/dm^^ [GeV-4] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0,85) 

2.00 X 10-3 

1.9 

7.6 

7.9 

[85,135) 

5.79 X 10-3 

1.7 

3.9 

4.3 

[135,190) 

4.44 X 10-3 

2.0 

3.0 

3.6 

[190,255) 

2.37 X 10-3 

2.5 

4.6 

5.3 

[255,325) 

1.13 X 10-3 

3.7 

7.6 

8.5 

[325,415) 

4.50 X 10-4 

4.9 

8.6 

10.0 

[415,505) 

1.57 X 10-4 

8.8 

10.9 

14.0 

[505,630) 

0.47 X 10-4 

12.6 

28.4 

31.1 

[630,800) 

0.12 X 10-4 

21.9 

15.4 

26.8 
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A Values of the normalized differential cross sections 


Table A.3: Normalized differential ft cross section in the dilepton channels as a function of the 
charged lepton transverse momentum (pj) and pseudorapidity {?](). The superscript 'i' refers 
to both and . The results are presented at particle level in the fiducial phase space. The 
statistical and systematic uncertainties are added in quadrature to yield the total uncertainty 


Pj bin range [GeV] 

l/crd<j/dpj [GeV 

Stat. [%] 

Syst. [%] 

Total [%] 

[20,40) 

1.93 X 10-2 

0.5 

1.6 

1.7 

[40,70) 

1.25 X 10-2 

0.5 

0.7 

0.9 

[70,120) 

3.79 X 10-3 

0.8 

1.8 

2.0 

[120,180) 

6.51 X 10-4 

2.0 

7.2 

7.4 

[180,400) 

3.89 X 10-3 

4.6 

12.0 

12.9 

bin range 

l/crd(7/dr]£ 

Stat. [%] 

Syst. [%] 

Total [%] 

[-2.4,-2.1) 

6.04 X 10-2 

3.0 

2.7 

4.1 

[-2.1,-1.8) 

9.49 X 10-2 

2.2 

2.9 

3.7 

[-1.8,-1.5) 

1.40 X 10-4 

1.9 

1.3 

2.3 

[-1.5,-1.2) 

1.98 X 10-4 

1.5 

1.5 

2.2 

[-1.2,-0.9) 

2.43 X 10-4 

1.4 

0.6 

1.5 

[-0.9,-0.6) 

2.86 X 10-4 

1.2 

2.3 

2.6 

[-0.6,-0.3) 

3.22 X 10-4 

1.2 

1.0 

1.5 

[-0.3,0.0) 

3.36 X 10-4 

1.2 

0.7 

1.4 

[0.0,0.3) 

3.17 X 10-4 

1.2 

1.2 

1.7 

[0.3,0.6) 

3.23 X 10-4 

1.2 

1.2 

1.7 

[0.6,0.9) 

2.90 X 10-4 

1.3 

1.1 

1.7 

[0.9,1.2) 

2.43 X 10-4 

1.4 

1.0 

1.7 

[1.2,1.5) 

1.97 X 10-4 

1.6 

1.2 

2.0 

[1.5,1.8) 

1.30 X 10-4 

2.1 

2.6 

3.4 

[1.8,2.1) 

9.70 X 10-2 

2.3 

3.1 

3.8 

[2.1,2.4) 

5.38 X 10-2 

3.5 

5.2 

6.3 
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Table A.4: Normalized differential ft cross section in the dilepton channels as a function of the 
transverse momentum (pj) and the invariant mass of the dilepton pair. The superscript 
't refers to both £+ and i~. The results are presented at particle level in the fiducial phase 
space. The statistical and systematic uncertainties are added in quadrature to yield the total 
uncertainty. 


Pj bin range [GeV] 

l/ada/dpj [GeV 

stat. [%] 

Syst. [%] 

Total [%] 

[0,10) 

1.84 X 10-3 

3.7 

3.2 

4.9 

[10,20) 

4.82 X 10-3 

2.0 

2.8 

3.5 

[20,40) 

7.74 X 10-3 

1.2 

2.8 

3.0 

[40,60) 

1.10 X 10-2 

1.0 

2.2 

2.5 

[60,100) 

9.70 X 10-3 

0.7 

0.9 

1.2 

[100,150) 

2.87 X 10-3 

1.3 

5.0 

5.2 

[150,400) 

1.03 X 10-4 

3.4 

7.0 

7.8 

m££ bin range [GeV] 

\l(jdcrldm££ [GeV-^] 

Stat. [%] 

Syst. [%] 

Total [%] 

[20,30) 

3.74 X 10-3 

2.6 

2.3 

3.5 

[30,50) 

5.28 X 10-3 

1.5 

2.1 

2.6 

[50,76) 

7.61 X 10-3 

1.0 

1.6 

1.9 

[76,106) 

7.19 X 10-3 

1.0 

2.1 

2.4 

[106,130) 

5.41 X 10-3 

1.3 

1.4 

1.9 

[130,170) 

3.30 X 10-3 

1.4 

1.7 

2.2 

[170,260) 

1.23 X 10-3 

1.5 

3.4 

3.7 

[260,400) 

2.29 X 10-4 

2.7 

6.9 

7.4 
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A Values of the normalized differential cross sections 


Table A.5: Normalized differential ft cross section in the dilepton channels as a function of 
the b jet transverse momentum (pj) and pseudorapidity (rjb), and the transverse momentum 
ip^) and the invariant mass of the bb system. The superscript 'h' refers to both b and b 
jets. The results are presented at particle level in the fiducial phase space. The statistical and 
systematic uncertainties are added in quadrature to yield the total uncertainty. 


Pj bin range [GeV] 

l/ada/dpj [GeV ^] 

Stat. [%] 

Syst. [%] 

Total [%] 

[30,50) 

1.16 X 

1.2 

7.7 

7.8 

[50,80) 

1.23 X 

1.2 

3.8 

4.0 

[80,130) 

5.94 X 10-3 

1.4 

3.4 

3.7 

[130,210) 

1.11 X 10-3 

2.2 

5.1 

5.5 

[210,400) 

6.88 X 10-3 

6.8 

13.5 

15.1 

t]\j bin range 

l/ad(7/dr]\, 

Stat. [%] 

Syst. [%] 

Total [%] 

[-2.4,-1.5) 

1.08 X 10-3 

2.3 

6.9 

7.3 

[-1.5,-1.0) 

2.16 X 10-3 

2.0 

1.8 

2.7 

[-1.0,-0.5) 

2.74 X 10-3 

1.9 

3.2 

3.7 

[-0.5,0.0) 

3.01 X 10-3 

1.9 

3.2 

3.7 

[0.0,0.5) 

3.20 X 10-3 

1.8 

2.6 

3.2 

[0.5,1.0) 

2.78 X 10-3 

1.9 

2.6 

3.2 

[1.0,1.5) 

2.19 X 10-3 

2.0 

1.9 

2.7 

[1.5,2.4) 

1.09 X 10-3 

2.4 

5.8 

6.2 

p^ bin range [GeV] 

l/ada/dp^ [GeV-3] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0,30) 

3.49 X 10-3 

2.3 

3.6 

4.3 

[30,60) 

6.50 X 10-3 

1.5 

2.9 

3.3 

[60,100) 

8.07 X 10-3 

1.2 

2.1 

2.4 

[100,180) 

4.27 X 10-3 

1.0 

3.2 

3.4 

[180,400) 

1.54 X 10-3 

4.4 

9.8 

10.7 

bin range [GeV] 

\l(j dcrldm^^ [GeV-3] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0,60) 

8.48 X 10-3 

1.9 

3.9 

4.3 

[60,120) 

4.57 X 10-3 

0.9 

2.7 

2.8 

[120,240) 

3.97 X 10-3 

0.6 

1.0 

1.1 

[240,600) 

5.35 X 10-3 

1.3 

4.6 

4.8 
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Table A.6: Normalized differential tf cross section in the £+jets channels as a function of top 
quark or antiquark observables: the transverse momentum (pj) and the transverse momentum 
in the tf rest frame (Pj) of the top quarks or antiquarks. The horizontally-corrected bin centres 
according to the MadGraph+pythia6 prediction (cf. Section |6.2| | are also provided. The 
results are presented at parton level in the full phase space. The statistical and systematic 
uncertainties are added in quadrature to yield the total uncertainty. 


Pj bin range [GeV] 

Bin centre [GeV] 

l/crda/dpj [GeV ^] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0,60) 

26.25 

4.14 X 10-3 

1.2 

3.6 

3.8 

[60,100) 

88.75 

6.69 X 10-3 

1.3 

1.7 

2.1 

[100,150) 

126.25 

4.96 X 10-3 

1.1 

3.0 

3.2 

[150,200) 

173.75 

2.66 X 10-3 

1.3 

3.5 

3.7 

[200,260) 

228.75 

1.06 X 10-3 

1.6 

3.2 

3.6 

[260,320) 

286.25 

3.99 X 10-^ 

2.2 

5.6 

6.0 

[320,400) 

356.25 

1.30 X 10-4 

2.8 

7.6 

8.1 

[400,500) 

446.25 

0.37 X 10-4 

5.5 

9.5 

10.9 

Pj bin range [GeV] 

Bin centre [GeV] 

l/crdcr/dpj [GeV-^] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0,60) 

26.25 

4.44 X 10-3 

1.6 

4.4 

4.7 

[60,100) 

63.75 

7.03 X 10-3 

1.7 

1.4 

2.2 

[100,150) 

126.25 

4.93 X 10-3 

1.5 

2.9 

3.2 

[150,200) 

173.75 

2.44 X 10-3 

1.9 

3.2 

3.7 

[200,260) 

226.25 

9.00 X 10-4 

2.4 

3.8 

4.4 

[260,320) 

286.25 

3.21 X 10-4 

3.3 

4.2 

5.4 

[320,400) 

356.25 

0.94 X 10-4 

4.5 

6.3 

7.8 

[400,500) 

443.75 

0.25 X 10-4 

9.6 

9.3 

13.4 
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A Values of the normalized differential cross sections 


Table A.7: Normalized differential tf cross section in the £+jets channels as a function of top 
quark or antiquark observables: the rapidity (i/t) of the top quarks or antiquarks, and the differ¬ 
ence in the azimuthal angle between the top quark and antiquark (A^(t,f)). The horizontally- 
corrected bin centres according to the MadGraph+pythia6 prediction (cf. Section |6.2| are 
also provided. The results are presented at parton level in the full phase space. The statistical 
and systematic uncertainties are added in quadrature to yield the total uncertainty. 


1/t bin range 

Bin centre 

1 l(T dcr/dl/t 

Stat. [%] 

Syst. [%] 

Total [%] 

[-2.5,-1.6) 

-2.01 

7.36 X 10-2 

2.5 

3.8 

4.5 

[-1.6,-1.2) 

-1.39 

1.75 X 10-1 

1.5 

2.4 

2.9 

[-1.2,-0.8) 

-1.01 

2.61 X 10-1 

1.4 

3.0 

3.3 

[-0.8,-0.4) 

-0.61 

3.00 X 10-1 

1.3 

1.3 

1.8 

[-0.4,0.0) 

-0.24 

3.33 X 10-1 

1.2 

1.5 

1.9 

[0.0,0.4) 

0.24 

3.31 X 10-1 

1.2 

1.6 

2.0 

[0.4,0.8) 

0.61 

3.00 X 10-1 

1.3 

1.8 

2.2 

[0.8,1.2) 

1.01 

2.47 X 10-1 

1.4 

3.0 

3.3 

[1.2,1.6) 

1.41 

1.88 X 10-1 

1.4 

1.6 

2.1 

[1.6,2.5) 

2.01 

7.77 X 10-2 

2.4 

3.4 

4.1 

A<^(t,f) bin range [rad] 

Bin centre [rad] 

l/crdo-/dA(p{t,l) [rad-i] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0.00,2.00) 

1.26 

6.83 X 10-2 

2.2 

7.0 

7.3 

[2.00,2.75) 

2.44 

3.22 X 10-1 

1.2 

3.8 

4.0 

[2.75,3.00) 

2.89 

1.13 

1.4 

3.2 

3.5 

[3.00,3.15) 

3.14 

2.27 

1.3 

6.6 

6.8 
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Table A.8: Normalized differential ft cross section in the ^+jets channels as a function of the px 
of the leading (pj) and trailing (pj) top quarks or antiquarks. The horizontally-corrected bin 
centres according to the MadGraph+pythia6 prediction (cf. Section [6^ are also provided. 
The results are presented at parton level in the full phase space. The statistical and systematic 
uncertainties are added in quadrature to yield the total uncertainty. 


Pj bin range [GeV] 

Bin centre [GeV] 

l/crdcr/dpj [GeV ^] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0,60) 

31.25 

2.61 X 10-^ 

2.2 

5.0 

5.4 

[60,100) 

76.25 

6.39 X 10-3 

1.5 

2.3 

2.8 

[100,150) 

126.25 

5.56 X 10-3 

1.3 

3.3 

3.6 

[150,200) 

173.75 

3.36 X 10-3 

1.6 

3.6 

3.9 

[200,260) 

228.75 

1.43 X 10-3 

1.9 

3.4 

3.9 

[260,320) 

286.25 

5.56 X 10-4 

2.5 

6.3 

6.8 

[320,400) 

356.25 

1.87 X 10-4 

3.1 

8.2 

8.8 

[400,500) 

446.25 

0.56 X 10-4 

6.0 

10.5 

12.1 

Pj bin range [GeV] 

Bin centre [GeV] 

l/crdcr/dpj [GeV^^] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0,60) 

23.75 

5.64 X 10-3 

1.4 

3.2 

3.5 

[60,100) 

83.75 

6.98 X 10-3 

1.9 

1.8 

2.6 

[100,150) 

123.75 

4.39 X 10-3 

1.7 

2.6 

3.1 

[150,200) 

173.75 

1.97 X 10-3 

2.2 

3.2 

3.9 

[200,260) 

228.75 

6.93 X 10-4 

2.8 

3.3 

4.3 

[260,320) 

286.25 

2.43 X 10-4 

4.0 

4.5 

6.0 

[320,400) 

356.25 

0.75 X 10-4 

5.3 

7.3 

9.0 

[400,500) 

443.75 

0.19 X 10-4 

11.4 

11.1 

16.0 
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A Values of the normalized differential cross sections 


Table A.9: Normalized differential tt cross section in the £+jets channels as a function of top 
quark pair observables: the transverse momentum (pj), the rapidity and the invariant 
mass im^l) of the tf system. The horizontally-corrected bin centres according to the Mad- 
Graph+pythia6 prediction (cf. Section |6.2| | are also provided. The results are presented at 
parton level in the full phase space. The statistical and systematic uncertainties are added in 
quadrature to yield the total uncertainty. 


Pj bin range [GeV] 

Bin centre [GeV] 

l/ada/dpj [GeV ^] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0,20) 

4.50 

1.32 X 

2.3 

7.9 

8.2 

[20,45) 

32.50 

1.18 X 10-2 

2.1 

3.4 

3.9 

[45,75) 

58.50 

6.40 X 10-3 

1.8 

5.5 

5.8 

[75,120) 

95.50 

2.84 X 10-3 

2.3 

6.6 

7.0 

[120,190) 

152.50 

1.07 X 10-3 

2.6 

5.7 

6.2 

[190,300) 

237.50 

3.06 X 10-4 

3.8 

12.2 

12.8 

J/tf bin range 

Bin centre 

1 1(7 da/dy^i 

Stat. [%] 

Syst. [%] 

Total [%] 

[-2.5,-1.3) 

-1.82 

6.07 X 10-2 

3.5 

2.9 

4.5 

[-1.3,-0.9) 

-1.09 

2.20 X 10-3 

1.8 

2.3 

3.0 

[-0.9,-0.6) 

-0.76 

3.27 X 10-3 

1.8 

2.9 

3.4 

[-0.6,-0.3) 

-0.46 

3.73 X 10-3 

1.7 

1.6 

2.3 

[-0.3,0.0) 

-0.14 

4.27 X 10-3 

1.6 

1.6 

2.2 

[0.0,0.3) 

0.14 

4.13 X 10-3 

1.6 

1.9 

2.4 

[0.3,0.6) 

0.47 

3.74 X 10-3 

1.7 

1.7 

2.4 

[0.6,0.9) 

0.76 

3.17 X 10-3 

1.8 

1.7 

2.5 

[0.9,1.3) 

1.09 

2.30 X 10-3 

1.8 

2.2 

2.8 

[1.3,2.5) 

1.82 

6.41 X 10-2 

3.3 

3.8 

5.1 

mjj bin range [GeV] 

Bin centre [GeV] 

l/crdcr/dwJtt [GeV-3] 

Stat. [%] 

Syst. [%] 

Total [%] 

[345,400) 

362.5 0 

4.69 X 10-3 

2.1 

7.1 

7.5 

[400,470) 

435.50 

4.30 X 10-3 

2.1 

2.9 

3.6 

[470,550) 

508.50 

2.67 X 10-3 

2.0 

6.1 

6.4 

[550,650) 

595.50 

1.17 X 10-3 

2.3 

7.3 

7.7 

[650,800) 

717.50 

4.66 X 10-4 

2.5 

4.2 

4.9 

[800,1100) 

927.50 

1.14 X 10-4 

3.0 

9.5 

10.0 

[1100,1600) 

1328.50 

0.11 X 10-4 

8.0 

9.8 

12.7 
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Table A.10: Normalized differential ft cross section in the dilepton channels as a function of top 
quark or antiquark observables: the transverse momentum (pj), the transverse momentum in 
the tt rest frame (pj), and the rapidity (i/t) of the top quarks or antiquarks, and the difference in 
the azimuthal angle between the top quark and antiquark The horizontally-corrected 

bin centres according to the MadGraph+pythia6 prediction (cf. Section |6.2| are also pro¬ 
vided. The results are presented at parton level in the full phase space. The statistical and 
systematic uncertainties are added in quadrature to yield the total uncertainty. 


Pj bin range [GeV] 

Bin centre [GeV] 

l/ada/dp\ [GeV ^] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0,65) 

28.75 

4.33 X 10-3 

1.1 

2.3 

2.5 

[65,125) 

101.25 

6.40 X 10-3 

0.9 

1.4 

1.7 

[125,200) 

161.25 

3.08 X 10-3 

1.0 

2.1 

2.3 

[200,290) 

238.75 

8.62 X 10-4 

1.3 

3.6 

3.8 

[290,400) 

336.25 

1.88 X 10-4 

2.7 

7.4 

7.9 

Pj bin range [GeV] 

Bin centre [GeV] 

1/ada/dpj [GeV-^] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0,60) 

26.25 

4.45 X 10-3 

1.5 

2.7 

3.1 

[60,115) 

93.75 

6.89 X 10-3 

1.4 

2.1 

2.5 

[115,190) 

151.25 

3.41 X 10-3 

1.3 

2.2 

2.6 

[190,275) 

226.25 

8.78 X 10-4 

2.2 

5.5 

5.9 

[275,380) 

318.75 

1.87 X 10-4 

2.6 

6.2 

6.7 

[380,500) 

428.75 

2.91 X 10-3 

7.5 

9.1 

11.8 

1/t bin range 

Bin centre 

l/adcr/dyt 

Stat. [%] 

Syst. [%] 

Total [%] 

[-2.5,-1.6) 

-2.01 

7.63 X 10-2 

2.0 

5.3 

5.6 

[-1.6,-1.0) 

-1.31 

1.97 X 10-4 

1.0 

1.5 

1.8 

[-1.0,-0.5) 

-0.76 

2.82 X 10-4 

1.1 

1.5 

1.9 

[-0.5,0.0) 

-0.29 

3.35 X 10-4 

1.0 

1.5 

1.8 

[0.0,0.5) 

0.29 

3.28 X 10-4 

1.0 

1.9 

2.2 

[0.5,1.0) 

0.76 

2.85 X 10-4 

1.1 

1.1 

1.5 

[1.0,1.6) 

1.31 

2.03 X 10-4 

1.1 

1.1 

1.5 

[1.6,2.5) 

2.01 

7.40 X 10-2 

2.0 

5.3 

5.7 

A<^(t,f) bin range [rad] 

Bin centre [rad] 

l/c7'd(7/dA<^(t,f) [rad-4] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0,1.89) 

1.19 

6.40 X 10-2 

2.0 

7.3 

7.6 

[1.89,2.77) 

2.44 

2.96 X 10-4 

1.0 

5.6 

5.7 

[2.77,3.04) 

2.94 

1.24 

1.1 

3.2 

3.3 

[3.04,3.15) 

3.09 

2.58 

1.5 

7.8 

8.0 
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A Values of the normalized differential cross sections 


Table A.ll: Normalized differential ft cross section in the dilepton channels as a function of the 
Pt of the leading (pj) and trailing (pj) top quarks or antiquarks. The results are presented at 
parton level in the full phase space. The statistical and systematic uncertainties are added in 
quadrature to yield the total uncertainty. 


Pj bin range [GeV] 

Bin centre [GeV] 

l/ada/dpj [GeV ^] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0,75) 

36.25 

3.31 X 10-3 

1.5 

3.7 

4.0 

[75,130) 

111.25 

6.43 X 10-3 

1.1 

1.8 

2.1 

[130,200) 

163.75 

3.67 X 10-3 

1.2 

2.9 

3.2 

[200,290) 

241.25 

1.17 X 10-3 

1.5 

3.9 

4.2 

[290,400) 

336.25 

2.61 X 10-4 

3.2 

6.9 

7.6 

Pj bin range [GeV] 

Bin centre [GeV] 

l/ada/dpj [GeV-^] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0,55) 

23.7 5 

5.38 X 10-3 

1.7 

2.2 

2.8 

[55,120) 

91.25 

6.74 X 10-3 

1.4 

1.6 

2.1 

[120,200) 

156.25 

2.50 X 10-3 

1.7 

2.1 

2.7 

[200,290) 

238.75 

5.58 X 10-4 

2.4 

4.6 

5.2 

[290,400) 

338.75 

1.14 X 10-4 

5.1 

10.2 

11.4 



47 


Table A.12: Normalized differential ft cross section in the dilepton channels as a function of 
top quark pair observables: the transverse momentum (pj), the rapidity (y^^) and the invariant 
mass (m^j) of the tt system. The results are presented at parton level in the full phase space. The 
statistical and systematic uncertainties are added in quadrature to yield the total uncertainty. 


Pj bin range [GeV] 

Bin centre [GeV] 

l/ada/dp^ [GeV ^] 

Stat. [%] 

Syst. [%] 

Total [%] 

[0,30) 

4.50 

1.43 X 

0.8 

6.1 

6.1 

[30,80) 

51.50 

6.90 X 10-3 

1.0 

4.7 

4.8 

[80,170) 

118.50 

1.91 X 10-3 

1.1 

5.4 

5.5 

[170,300) 

223.50 

3.47 X 10-4 

2.1 

4.7 

5.1 

1/tt bin range 

Bin centre 

1 Id dcr/di/tt 

Stat. [%] 

Syst. [%] 

Total [%] 

[-2.5,-1.5) 

-1.93 

4.71 X 10-2 

4.0 

6.6 

7.7 

[-1.5,-1.0) 

-1.24 

1.82 X 10-4 

1.8 

1.5 

2.3 

[-1.0,-0.5) 

-0.76 

3.09 X 10-4 

1.4 

1.5 

2.1 

[-0.5,0.0) 

-0.31 

4.18 X 10-4 

1.2 

1.7 

2.1 

[0.0,0.5) 

0.29 

4.09 X 10-4 

1.2 

1.2 

1.7 

[0.5,1.0) 

0.76 

3.15 X 10-4 

1.5 

1.3 

2.0 

[1.0,1.5) 

1.24 

1.79 X 10-4 

1.8 

1.8 

2.6 

[1.5,2.5) 

1.93 

4.59 X 10-2 

4.0 

6.0 

7.2 

bin range [GeV] 

Bin centre [GeV] 

l/trdn/dmjj [GeV-4] 

Stat. [%] 

Syst. [%] 

Total [%] 

[340,380) 

354.50 

4.14 X 10-3 

3.0 

8.6 

9.1 

[380,470) 

428.50 

4.50 X 10-3 

1.7 

5.3 

5.6 

[470,620) 

537.50 

1.95 X 10-3 

1.8 

2.9 

3.4 

[620,820) 

705.50 

5.25 X 10-4 

2.8 

3.2 

4.2 

[820,1100) 

940.50 

1.00 X 10-4 

3.7 

7.3 

8.2 

[1100,1600) 

1328.50 

7.28 X 10-6 

14.4 

28.2 

31.6 
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A Values of the normalized differential cross sections 
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